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Abstract 

We study evolution of entanglement of two two-level atoms placed inside a multilayered mi- 
crosphere. We show that due to inhomogeneity of the field modes this entanglement essentially 
depends on the atomic positions (asymmetrical entanglement) and also on the detuning between the 
atomic transitions and field frequencies. The robust and complete entanglement can be achieved 
even in the resonant case when the atoms have different effective coupling constants, and it can 
be extended in time if the detuning is large enough. We study analytically the lossless case and 
estimate numerically the effect of dissipative processes. 
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I. INTRODUCTION 



Recently, essential progress in fabrication and determination of optical properties of dif- 
ferent kinds of microcavities with sizes about 0.1 — 20fim which contain semiconductor 
nanoclusters or quantum dots (QDs) has been achieved [see and references therein]. 

When semiconductor QDs are embedded in a spherical microcavity, the QD luminescence 
can be coupled with eigen modes of the electromagnetic field of the microcavity and a lower 
ihreshold of stimulated emission (or lasing modes) of QDs can be achieved. In recent papers 

inn 

3],|3|,|^ a coupling between the optical emission of embedded CdSe^Si^x , QDs and spher- 
ical cavity modes was studied and a strong whispering gallery mode (WGM) resonance with 
high Q factors is registered in the photoluminescence spectra. Recently quantum-confined 
semiconductor nanorods were used as highly polarized nanoemitters for active control of the 
polarization state of microcavity photons. 

Until now the modes with small numbers of spherical harmonics (SNM) are essentially 
less well-studied compared with whispering gallery modes (WGM) due to their rather low 
Q-factor caused by significant radiating losses. A possibility of strongly increasing the Q- 



facto, of tKe „.c..pKe.e wa. ,opo.ed . several pape., see e., fl, ^ fl. TKe .a. .ea 

consists in coating a microsphere by alternative layers of a spherical stack, which results in 
an increase of the Q factor up to values comparable for WGM, i.e. 10^ — 10^. 

In a system with small mean photon number two spatially separated atoms in a cavit 
become entangled at some time moments as a result of sharing the re-radiated photons 
3. 0. Q. Q. one of these sche.es has been .called us,ng Rvdbcg ato,„s coupled 
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one by one to a high Q microwave superconducting microcavity [1^. In inhomogeneous 
structures, like multilayered microspheres, the quantized field properties are quite different 
from the unbounded case because of a non-uniformity of the cavity field, which becomes 
important for the entanglement dynamics. In spite of numerous experimental obstacles, 
mainly related with the decoherence problem, it seems very natural to entangle atoms placed 
in high-Q cavities (like microspheres) via interaction with modes of the cavity quantized field. 
A simple scheme for the generation of two-atom maximally entangled states via dispersive 
interaction was proposed in A number of papers report studies of the evolution of the 
entanglement in an atomic subsystem resonantly interacting with a single mode of the cavity 



onantiv inters 



field (two-atom Tavis-Cummings model) [la], [131, [iSl- A robust generation of many-particle 
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entanglement in various configurations has been discussed in Q], [2^, j^|. Experimentally 
a robust entanglement was recently studied inj^ where an entanglement lasting for more 
than 20s was observed in a system of two trapped Ca+ ions. Authors jssi have shown that 
the degree of entanglement between the two atoms strongly depends on the mean photon 
number and the strength of two-photon correlations. 



In 



2J| a scheme for entangling two-level atoms located close to the surface of a dielectric 



microsphere and atoms resonantly interacting with the field was considered. It was shown 
that in the particular case of two atoms located at diametrically opposite positions a perfect 
entanglement cannot be achieved even in the strong-coupling regime. 

In this paper we study two-atom entanglement interacting with field modes inside a 
microsphere covered with spherical dielectric alternating layers (coated microsphere). We 
are mainly interested: (i) m the frequency range of the high reflectivity field in A/4 -stack; 
(a) in the case when identical atoms are located asymmetrically inside the microsphere 
(i.e. the system is not symmetric with respect to a permutation of initially excited and 
non-excited atoms), so that the field inhomogeneity leads to different effective atom-field 
coupling constants ; and (in) the atomic transitions can be both resonant and well detuned 
from the field peak frequency. 

The paper is organized as follows. In Section II we discuss basic equations for two atoms 
placed into a coated microsphere and the solution for this case. In Section III we present an 
analytical solution for probability amplitudes and apply it to studying the atomic concur- 
rence. In Section IV we present a numerical study of the concurrence (tangle) dynamics. In 
the last Section, we discuss and summarize our conclusions. 



II. BASIC EQUATIONS 

Consider two identical two-level atoms coupled to a quantized electromagnetic cavity field 
in a coated microsphere (Fig^. 

Let us assume that the atoms are sufficiently far from each other, so that the interatomic 
Coulomb interaction can be ignored. In this case, the electric dipole and rotating wave 
approximations can be applied and the Hamiltonian for the atom-cavity system [h = 1) is 
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FIG. 1: Geometry of coated microsphere with two atoms. 



given by 




H = Ho + H,, (1) 
J Jo ^.^1^2 ^ 

j 

where Uj is the atomic transition frequency (cui = U2 = ujat) , Sz,±j, j = 1,2 are the 
atomic operators corresponding to the j-th atom and obeying standard su{2) commutation 
relations, [si,^^] = ±s±, = 2sz, dj are atomic dipoles. Here f{r,uj) and f'^{r,u) are 

bosonic operators which play the role of the fundamental variables of the electromagnetic 
field and the medium, including a reservoir necessarily associated with losses in the medium. 



The electric-field operator is expressed in terms of f{r,uj) as 2^, 

^(+)(r) =^./A Hdu^ [ d3r'v/i^(^G(r,r',.;)/(r', 

V TTEo Jo J 



(2) 



with 



Mr, uj), fUr', Lo') = 5,,5(r - r')6{uj - u'), 



(3) 



/i(r,'^),/i(r',^') 



/J(r,c.),/t(r>') 



where G{ r, r', u) is the classical Green tensor satisfying the equation 



e( r, Lu) — V X Vx 



G( r, r\uj) = -S{r- r') 



(4) 



together with the boundary condition at infinity [S{r) is the dyadic ^-function]. Here e{r, u) = 
eji{T,uj) + iei{r',u) is the complex dielectric permittivity. We look for the solution of the 
Schrodinger equation with the Hamiltonian in a single excitation manifold in the form 



|vl>(t)) = C^{t) |0) \e^g2) + C^it) |0) \g^e2) + 

//•OO 
d^r duj[Csiir,uj,t)\{hir,uj)})], 



(5) 



where le^) {\gk)) denotes the excited (ground) atomic state of k-th atom. Correspondingly, 
\{li{r,u)}) = f- {r,uj)\{0}) is a single photon Fock state and, |{0}) is the vacuum state of 
the rest of the system. Note that this state is not a photonic state in g eneral, but a state of 



24|. 



the macroscopic medium dressed by the electromagnetic field 27|.|28|. 2^ 

For simplicity we study the frequency range close to the microsphere resonance with the 
frequency ujf, when the Green function can be written as 



^(r, r', u) = G{r, r', u) ■ 6{uj — ^f) 



(6) 



The effect of broadening of such a line due to dissipation is studied in Sec. IV numerically. 
For further references we show in FigOthe typical frequency spectrum of Green's function, 
calculated numerically. Let us assume that the atomic dipoles are parallel to the surface of 
the microsphere (similar to the situation considered in 5[), so only tangential components 
of Green's tensor (e.g. Gyjy?) give a contribution. 

Projecting |\E'(t)) in © onto |0) \eigk) and |{lj(r, c<j)}) \g1g2) states, we obtain the follow- 
ing equations for the probability amplitudes 



(7) 



B{r,ujf,t) = iAujB{r,ujf, t) — iG(r, Oi, a;/) Ci{t) — 'iG(r, 02, cu/) C2(t), 
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FIG. 2: (a), (b) and (c). Frequency spectrum of imaginary parts of tangential component of the 
dyadic Green's function lm{G^^{r, r', /)),/ = ci;/27r for 7-layered system (microsphere coated with 
5 alternating A/4 layers), with atomic positions ai = 0.9/im and a2 = l.l^m. Refraction indexes of 
the layers are n4 = 1.5 + i2-10~^ (glass, bottom microsphere, Ifim), = 3.58 + 110"^ {Si, 0.12/im), 
712 = 1-46 + i3 ■ 10~^ {Si02, 0.3/um) and ni = 1 (surrounding space), (a) Im(G^^(ai, ai, /)); (b) 
Im(G(p^(ai, 02, /)); and (c) Im(G^^(a2, 02, /)); (d) radial dependence of Im(G<^^(r, 02, //)), where 
the atom is placed in 02 = 0.9^m and the field's peak frequency is // = 24:1.7THz. Dash line in 
(d) shows the refraction indexes of the spherical stack structure. 

where r is coordinate vector, 01,2 are the positions of the atoms in the microsphere, Au; = 
Uf - Uat, and 

B{v,ujf,t) = di j d^r'-aGjfc(a,r',a;/)C3fc(r',a;/,t), (8) 

G(a, r', uj) = ndidk Im(Gifc(a, r', lu)), (9) 

where a = i^J e i{y , uj f) / t^Squj'^j / and k = uP'^j (^t^Ei^. Eliminating B{r,ujf,t) from ((Tj) we 
obtain after minor algebra closed equations for Ci^2{t) in matrix form as follows 

d'^o . da ^ , . 

*T-^"i + A.q=0. (10) 
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where 



q 







, A = 


C2 





Im Gki{r,rl,u) 



G(l,l) G(l,2) 
[g(2,1) G(_2,2) 
To derive (jlOj) the identity 

/ d^s (tu^/c^) £1(8, a;)G'fcm(r, s, co')G'^^(r', s, cu) was taken into account. From now on 
we adopt the convention of summation over repeated vector-component indices. 
The general solution of Eq. ljlOj) has the form 

4 

Cfc(t) = ^Cfc,e'"^*, fc = l,2, 



where the frequencies 



det [(— cu^ + coiAco") 5ki + Aki] = or 



(11) 



are solutions of the eigenvalue problem 



+ uj/\ujY + (-u;^ + uAuj) Tr{A} + det(A) = 0. (12) 
Following [2^ we rewrite the Green tensor for a multilayered microsphere as follows 

G(r, r', u) = G^(r, r', u)5fs + G^^^) (r, r', u), (13) 

where G^{r, r', u) represents the contribution of the direct waves from the radiation sources 
in an unbounded medium, / and s denote the layers where the field point and source point 
are located, Sfs is the Kronecker symbol, and the scattering Green tensor G^-^*-*(r, r', cu) 
describes the contribution of both multiple reflection and transmission. The Green tensor 
q(/s) general can be expanded as 



Qife)(^ r' cu) = — W V + l i^-^)- (2-Sr, )G(-^'^)fr r' 
^ ^^'^''^^ 4^ 2^ 2^ 2^ n{n + l) {n+m)\^ Oon^j W„„.lr, r , 



(14) 



where G^cnm (r, r/, lS) is a particular Green tensor, n is the spherical and m is the azimuth 



quantum numbers of a microsphere, fcj = ujrii/c, rii = ^^/ei(iu) is a refraction index. General 
recurrent formulas and particular representations of the Green tensor Gpnm{r,r/,uj) can be 
found in Ref.H- 

Generally analysis of the Green tensor p4j) requires intensive computation. In the sim- 
plest case when atoms are located at positions with the same value of the field ampli- 
tude, we have G{i,j) = G and det(A) = 0. This case is symmetrical with respect to 



permutation of the atoms and one can easily obtain the solution ui = 0, uj2 = Auj and 
ci;3,4 = {Alj/2) ± [(Aco'/2)^ + 2G]^''^. However, experimentally, a symmetric location of the 
atoms with respect to the center of the microsphere is difficult to achieve In the sim- 
plest nontrivial case we have to take into account the nonuniformity of the field. In this case 
we have G{1, 1) 7^ G{2, 2) , but the coefficients G{i, k) can be written as follows 

1) = xl G{2, 2) = xl 2) = G(2, 1) = X1X2, (15) 

so that the condition 



G(1,1)-G(2,2) = G(1,2)2 (16) 
is fulfilled. We have again det(A) = 0, but now u^^^ 

G{1, 1) + G(2,2), so that the atoms have different Rabi frequencies. Introducing the new 
variable C3 according to 

Cs - lAuCs = -I (xiCi + X2C2) , (17) 
the system can be reduced to a simple form 

Ci = -^XiCs, C2 = -^XsCa- (18) 

From direct calculations we have found that such a case is fulfilled for a spherical structure 
with a 7-layered system (microsphere coated with 5 alternating A/4 layers), see FigI21with 
i?2 = 1.97/im, Ri = Ifim. Two atoms having tangentially oriented dipoles d ± r are at 
positions ai = O.dfim and 02 = l.lfim correspondingly. In FigI2fd) we show the radial 
distribution of the Green tensor component Im(G<^^) for the first atom in the microsphere. 
We have found that for this structure Xi = 1-72, X2 = 0.608 and X2/X1 = 0.35. As Xi > 
X2 we can see that in this case the coupling constant is larger for the first atom or, in other 
words, the first atom is placed in a stronger field mode of the coated microsphere. 
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III. EFFECTIVE HAMILTONIAN DYNAMICS 



In the lossless case the Hamiltonian corresponding to the simplified situation in the single- 
mode regime described by Eq.® can be represented in the following form 

H = ujftt^a + LUat {s^i + 5^2) + Xi (as+i + h.c.) + X2 + h.c.) , (19) 

where the effective coupling constants Xi = i)^^"^ depend on the positions of atoms inside 
the microsphere. Because the coupling constant is larger for the atom placed in the region of 
a stronger field mode, the configuration is not symmetrical with respect to the permutation 
of the atoms. In this case the state vector is given by 

\^{t)) = C,{t) |0) \e,g2) + C2{t) |0) \g,e2) + Csit) |1) \gig2) , (20) 
where Ci^2{t) are solutions of the Eqs.(fT7 j) .(fTH j) : 



Ci{t) = -xAt) + \, C2it) = -X2r{t) + 1-X, Cs{t) = -zixxMexp{zAut/2)sm{nt), (21) 



with 



Xx = XiA + X2(l - x'^ = xl + xl 



*^?T sm(iZt) — cos{\ It) 



+ 1}, 



(22) 



and the initial conditions Ci(0) = A, 6*2(0) = 1 — A, (A = 1, 0) are considered (for the A = 
case the system evolves from the initial state \g1e2)) ■ In particular, the average photon 
number can be easily calculated using the solution (j2T| : 



The reduced atomic density matrix for the state (j20|l has the form 



p- = Trf{\^) {^\} 











\Cif \C1C2 

IC1C2I |C*2| 

















(23) 
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In the frame of the standard approach j30| we obtain from ()23|) the concurrence C{t) for 
two atom system as 

C{t) = 2\CiC2\. (24) 

In Fig|21 we show the dynamics of tangle (see Eq.(j2D)) for the cases A = (Fig|21^a)), 
and A = 1 (FigOl^b)). For this configuration the condition Xi > X2 is fulfilled, in FiglHl^a) 
the first atom being in the ground state is placed in a larger field strength. In this case 
the tangle has the form of well resolved periodical plateau. The inverse situation is 
shown in FiglSfb). In Fig|21the dashed line shows the average photon number (n). One can 
observe that the amplitude of rapid oscillations in the upper part of a long periodical tangle 
evolution (plateau) is essentially less in the case when an unexcited atom is placed in 
the region of stronger field (FiglH^a)). It is easy to see, that such oscillations have the Rabi 
frequency and are related to the instantaneous average number of photons (dashed line in 
Fig|2)) stored in the field. 

It is worth noting that in the case when (7(1,1) • G{2,2) 7^ G(l,2)^ the general for- 
mulae (fTUI) should be used to study the dynamics of concurrence C. The evolution of the 
concurrence for the general situation is shown in Figl^ (solid line). The particular situation 
described by the simplified model (PT|) is presented in the same Figl^ as the dashed line. 
From FigEJone can observe that if (fTBj) is fulfilled the amplitude of the fast oscillations is 
much less than in the general case, and therefore can be regarded as the optimal dynamics. 

For the far detuned case \Auj\ » x can easily obtain from (PT|) two well sepa- 
rated frequency components of Ci,2(^):the high frequency component exp(zAti;t) and the 
low frequency component exp{ig^t / Au) . The latter leads to a formation of well recognized 
plateaus, which do not exist in the resonant case. Further, we will use the initial conditions 
corresponding to A = 1, so that Xx = Xi- The concurrence for such a solution can be 
easily analyzed in parameter space XiJ X2 f*^^ ^^e case Auj = 0. In this case concurrence ()24|) 
is explicitly asymmetrical with respect to the atomic permutation {xi ^ X2) has the 
form 



Cixi,X2) = 2^'^' 



1 - ^ [1 - cos(fit)] 



[1 - cos(f]t)] . (25) 



^]2 

It is easy to see from that the surface C{xi,X2) is separated by circles (with radii 
(Xi + xl)^^'^ = 2kiT/t, k = 0,1, 2..) on which C{xi, X2) = (atoms are disentangled). The 
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FIG. 3: Dynamics of the two-atom tangle 



vs r = ujatt for Xi = 0.254 and X2 — 0-151 in a 



lossless case for (a) initial state \g\e2]\ (b) initial state \e\g2]- Due to large detuning IS.uj/ujat — 0.5 
the amplitude of mean photon number (n) oscillations is less in (a) case. One can see the well- 
recognized plateaus of tangle in case (a). In case (b) the plateaus is adding with strong oscillations 
due to large (n). 

detailed structure of the concurrence can be better understood rewriting Eq. ()25|) in the form 



where k = X2/X1 the quantity a = 1 — cos(fit) is in the range < a < 2. At a 
fixed value of a the concurrence C{k,a) assumes maximal values C{k,a) < 1 at ki^2 = 



Note that k2 exists for a > 1. In general 1 < A;i < ^ 2.41 and < A;2 < A;^ ^ 0.41. This 
means that in the resonant case, Au = 0, the two-atom system can be maximally entangled 
if X2/X1 = ! when the atoms have different field-atom coupling constants, X2 Xi- 
Nevertheless in the general case, when Auj ^ , the structure of the concurrence C{xi, X2) 
is more complicated. 

The structure C(xi,X2) calculated for Uatt = 27, A = (initially excited atom is 




(26) 



2-1 6a±2(9a2 -4a + 4) 



1/2 



, and C{k,a) = 1 only at a = 2 when k"^ = v^± 1 = 1/k^. 
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FIG. 4: Tangle vs r = coatt in general case for parameters ALo/tOat = 0.75 , and (a) G{1, 1) = 
0.01, G(l,2) = G(2,l) = 0.012 and G(2,2) = 0.04. Dash line shows tangle for G(l, 1) = 0.01, 
G(l,2) = G{2, 1) = 0.02 and G(2,2) = 0.04 when Eq.JTHl) is valid. In case (b) G(l, 1) ^ G{2,2). 



placed in a smaller field) and Au/ujat = 0.5. We observe that the surface C(xi, X2) is rather 
asymmetrical with respect to the line Xi = X2- the course of evolution for fixed X17X2 
the maximal values of C{xi, X2) move out from the origin of coordinates. Obviously on the 
edges where Xij X2 = 0^ ^he concurrence vanishes, C(xi, X2) = 0- the vicinity of maxima 
the concurrence C{xi,X2) is highly asymmetric. For Xi > X2 the maxima C{xi,X2) i'^ 
the left side are smoother and the hills are more pronounced. This means that the system 
remains in the region of strong entanglement for long periods if Xi > X2- However in general, 
the details of the surface C(xi,X2) essentially depend on Au and the form of the Green's 
function. 

IV. NUMERICAL STUDY 

It is worth noting, that in a real microsphere the field dissipation is caused by material 
losses and the radiation into surrounding space leads to line broadening (bandwidth). The 
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FIG. 5: The same as in FigEl but for loss case (7 = 2 • 10 ^). To see the details of long-time 
dynamics we calculate up to Tmax = 1000. 

analytical calculation of such a broadening requires an extensive knowledge of the microscop- 
ical local field, which in a multilayered microsphere case is itself a quite difficult problem. 
To estimate the influence of the dissipation on the concurrence dynamics we will use the 
following simplifled approach. Although we do not know the exact frequency dependence 
of the dissipative part on the refractive indices of the materials in a microsphere, it is 
possible to calculate the spectral width of the Green function peak. Thus, we can estimate 
the effect of the fleld's dissipation using the master equation technique in the framework 
of the Lingblad approach. In particular, the dissipation coefficients are calculated from the 
bandwidth of the Green function peak (see FigH]). Such a semi-analytical approach allows 
us to simulate numerically not only the evolution of the concurrence in a lossy environment, 
but also the dynamics of the average photon number. In this approach we replace the ex- 
act Hamiltonian by the simplifled Hamiltonian ^T^i and numerically solve the following 
master equation for the joint atom-fleld density operator p in a dissipative cavity at zero 
temperature: 
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i[H, p] + Lip, 



+ 



a ap — pa a 



(27) 



(28) 



where H is given by ()19p. and in ()27p we neglect the atomic spontaneous emission in the 
Rabi period time scale. Also we have used the detuning Au/uat = 0.5, and Xi = 0.254, 
X2 = 0.151, X2/X1 = 0.594. In FigElthe dynamics of the two-atom tangle for the lossy 
case is shown. It is clear from Figl^l that the plateaus of concurrence survive even in the 
presence of dissipation, although their amplitude is obviously lower than in the lossless case. 

V. CONCLUSION 

In conclusion, we have studied the dynamics of entanglement of spatially separated two- 
level atoms interacting with a radially nonuniform cavity field mode in a dielectric micro- 
sphere coated with an alternating stack. We found that due to the field inhomogeneity the 
atoms can be maximally entangled even in the resonant case. We have found that entangle- 
ment essentially depends on the atomic positions (asymmetrical entanglement) and also on 
the detuning between atoms and the field mode frequencies. The entanglement is consid- 
erably more stable with duration much longer than the period of Rabi oscillations (robust 
entanglement) when the unexcited atom is placed in a stronger field, while the detuning 
increases the duration of the entanglement period. The dissipation reduces the amplitude 
of the entanglement, however practically does not change the width of the zones of large 
entanglement. 
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